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ABSTRACT 

Astro-H wiW be able for the first time to map gas velocities and detect turbulence in galaxy clusters. One of the best 
targets for turbulence studies is the Coma cluster, due to its proximity, absence of a cool core, and lack of a central 
active galactic nucleus. To determine what constraints Astro-H will be able to place on the Coma velocity field, we 
construct simulated maps of the projected gas velocity and compute the second-order structure function, an analog of 
the velocity power spectrum. We vary the injection scale, dissipation scale, slope, and normalization of the turbulent 
power spectrum, and apply measurement errors and finite sampling to the velocity field. We find that even with sparse 
coverage of the cluster, Astro-H will be able to measure the Mach number and the injection scale of the turbulent power 
spectrum-the quantities determining the energy flux down the turbulent cascade and the diffusion rate for everything 
that is advected by the gas (metals, cosmic rays, etc.). Astro-H will not be sensitive to the dissipation scale or the 
slope of the power spectrum in its inertial range, unless they are outside physically motivated intervals. We give the 
expected confidence intervals for the injection scale and the normalization of the power spectrum for a number of 
possible pointing configurations, combining the structure function and velocity dispersion data. Importantly, we also 
determine that measurement errors on the line shift will bias the velocity structure function upward, and show how 
to correct this bias. 

Subject headings: galaxies: clusters: intracluster medium — techniques: spectroscopic — X-rays: 
galaxies: clusters — methods: numerical 


1. INTRODUCTION 

X-ray observatories have yielded a wealth of informa¬ 
tion about the thermodynamic and chemical properties 
of the intracluster medium (ICM) of galaxy clusters. 
One aspect of the ICM that has been up to now 
beyond the ability of present instruments to directly 
measure is its kinematics. Determining the kinematic 
properties of the ICM is essential for a complete picture 
of the physics of the cluster gas. Kinetic energy in 
the form of bulk motions and turbulence likely pro¬ 
vides non-negligible pressure support against gravity, 
biasing mass estimates based on the assumption of 
hydrostatic equilibri u m, as predi c ted b y simulations 
(lEvrard et al.l [T 996I iRasia et al. I2006I iN agai et ^ 
l2007t iPiffaretti fc Valdarninil 12008 iNelson et al.11201^ 
possibly explaining discrepancies bet ween hydro¬ 
static and weak lensing-deriyed masses (IZhang et al.l 
20IOI: [Mahdayj et al.l l2013t Ron der Linden et al.l 12014 


Applegate et al.l 1201411 . Dissipation of turbulent kinetic 


energy into heat and turbulent transport and mixing 
of hot gas m ay partially o ffset g a.s cooling in cluster 
cool cores (iFuiita et 12004 iDennis fc Chandrad 
20051 l^uHone et al.l 2010l : iBaneriee fc Sharmal 12014 


Zhurayleya et al.ll2014bll . The properties of gas motions 


place important constraints on the microphysics of the 
ICM, in particular its yiscosity (iFabian et all 1200.4 
iRoediger et al.l[2?)T4 iZuHone et al.ll2015h . Finally. ICM 
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turbulence is likely a key ingredient for the origin of 
non-thermal phenomena such as radio h alos and radio 
mini-halos (iBrunetti fc: LazarianI 120071: iDonnert et alJ 
IW^IZiiHouc et RlJl2ni,HL 

Though certain o bseryations, such as that o f cold f ronts 
(see the reyiew by iMarkevitch fc VikhlininI (l2007fl and 
the re cent simulation inyestigati ons by Roediger et al.l 
(|2013f) and iZuHone et al.l (|2015ll for details), indicate 
that the ICM may be somewhat yiscous, the Reynolds 
number is probably high enough to permit the deyel- 
opment of turbulence. Major and minor mergers with 
other clusters, actiye galactic nucleus (AGN) outbursts, 
and stirring of the ICM by the motion of cluster galax¬ 
ies and dark matter substructure are all possible driyers 
of turbulent gas motions. In the simplest picture, tur¬ 
bulent motions are driyen at a largest, “injection” scale, 
and cascade down to smaller and smaller scales until they 
are dissipated at the smallest or “dissipation” scale. 

In this picture, the distinction between “bulk flows” 
and “turbulence” is blurry-most merger-induced large- 
scale bulk flows simply represent turbulence at its injec¬ 
tion scale. A “bulk flow” in a merging cluster could not 
be considered part of the turbulent field if the moying gas 
parcel is still grayitationally bound to an infalling sub¬ 
cluster (i.e., it is still part of the “spoon” that stirs the 
turbulence); when that gas is stripped by ram pressure, 
it joins the turbulent yelocity field. Infalling subclusters 
are usually easy to identify in the X-ray images and, if 
necessary, can be masked for a study of the turbulent 
yelocities. 

Physical models of the plasma combined with dimen¬ 
sional arguments yield simple power-law forms for the 
turbulent power as a function of wayenumber, the most 
familiar of which is the Kolmogoroy spectrum. There¬ 
fore, obseryationally determining the relationship be¬ 
tween measured gas yelocities and the range of length 
scales in the cluster has the potential to reyeal the un- 
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derlying gas physics. In particular, determining the dis¬ 
sipation scale of the turbulent cascade would give the 
effective “viscosity” of the ICM, which is governed by 
complex dissipation processes in the magnetized, col¬ 
lisionless plasma. Knowledge of the injection scale of 
the cascade, in combination with the average veloc¬ 
ity dispersion (or Mach number) of the turbulent mo¬ 
tions, would give us the diffusion rate of everything 
that is advected by the gas (such as metals, cosmic 
rays, specific entropy. rRebnsco et al.l 1200^ iVazza et al.l 
I 2 OIOI : lEnfllin et al.M201lll . because the largest-scale ed¬ 
dies are those most responsible for stirring and mixing 
the ICM. This combination also determines the flux of 
energy that moves down the turbulent cascade and ends 
up dissipating into heat and accelerating cosmic ra ys 
(jPennis fc Ctiandranll2005l iBrunetti fc Lazarianll2007fl . 

Turbulence i n the ICM can be studied dir ectly by sev¬ 
eral methods. Ilnoeramov fc SnnvaevI (1200,111 pointed out 
that the finite number of the largest-scale turbulent ed¬ 
dies that fit in a cluster volume will result in devia¬ 
tions of the spectral line profiles from a simple Gaus¬ 
sian shape. This method requires well-resolved spectra 
with high statistical quality to detect subtle shape de¬ 
viations. Another method relies on mapping the pro¬ 
jected line of sight velocities over the face of a cluster 
to derive a velocity structure function, which has a one- 
to-one correspondence t o the turbulent power spectrum 
(j Zhuravleva et all 120121 hereafter Z12, and references 
therein). This method is less demanding to the statis¬ 
tical accuracy of the line profiles-only the line positions 
are used. We will consider this method below. 

So far. X-ray observatories have lacked the required 
spectral resolution to measure the line shifts and broad¬ 
ening resulting from turbulent motions in the ICM. The 
RGS grating on XMM-Newton can provide weak up¬ 


per limits on Do ppler broadening of spectral lines in 
cool- c ore clusters (ISanders et al.ll201lHSanders fc FabianI 
I 2 OI 3 I : iPinto et al.l 120151. and references therein). In¬ 
direct estimates of the ICM turbulent velocity can 
be obtaine d from measur e ments of resonant s c atter¬ 
ing (e.g.. iChurazov et ail 12004 IWerner et al.l I2009I 


Ide Plaa et M .1 I 2 OI 2 I: iZhuravleva et al.l 1201, III , pressure 
fluctuations (|Schuecker et al. [HoH o r surface bright¬ 
ness fl uctuations ( Churazov et al.ll2012l : IZhuravleva et al.l 


The Astro-H Mission, a joint JAXA/NASA endeavor, 
will be launched in late 2015 and will for the first time 
have the energy resolution to detect line shifts and widths 
resulting from bulk an d turbulent motions in the ICM 
(|Takahashi et al.l 1201211 . Astro-H will possess a Soft X- 
ray Spectrometer (SXS) micro-calorimeter with an en¬ 
ergy resolution of AA < 7 eV within the energy range 
E ^ 0.3— 12.0 keV, covering a 3’x3’ field. At the energy 
of the Fe-Kc line, E k, 6.7 keV, this enables the mea¬ 
surement of velocities at resolutions of tens of km s“^. 

A promising target for measuring the kinematic prop¬ 
erties of the ICM with Astro-H is the Coma cluster. 
Coma is nearby (^ 100 Mpc, 0 = 0.0231), possesses a flat 
core with a radius of Tc ^ 300 kpc (Rill’), and is likely 
to be turbulent due to ongoing merger activity. Further, 
Coma lacks a central AGN, which often generate ICM 
motions in their vicinity. Its lack of strong entropy gra¬ 
dients, unlike those in cool-core clusters, implies that the 
effects of stratification are not significant, and turbulence 


should develop in an isotropic, qualita tively simple way. 
We do note that ISanders et al.l (|2013[1 argued that the 
presence of large-scale filamentary structures they have 
detected in the X-ray image of Coma implies that tur¬ 
bulence is suppressed. However, such large-scale coher¬ 
ent structures are seen in hydrodynamic simulations in 
turbulent clusters fe.g. IVazza et al.ir2009ll . so a compari¬ 
son with simulations is needed to draw conclusions from 
those observations. 

In this work, we will simulate turbulent velocity fields 
in the Coma cluster, using a simple form for the under¬ 
lying power spectrum to test the ability of Astro-H to 
constrain the slope of the power spectrum and relevant 
length scales, such as the injection scale of the turbulent 
motions and the dissipation scale. From these simulated 
velocity fields, we will construct the second-order struc¬ 
ture function, a measure of the spatial structure of the 
velocity field that is equivalent to the power spectrum, 
and fit it with various models, in combination with the 
velocity dispersion. We will investigate a number of pos¬ 
sible spatial configurations for individual SXS pointings, 
in order to determine which are most useful for constrain¬ 
ing particular aspects of the velocity field. We assume a 
flat ACDM cosmology with h — 0.71 and Pm = 0.27. At 
the distance of Coma, an angular size of 1’ corresponds 
to a length scale of ~26 kpc. 

2. METHOD 

In what follows, we closely follow the work of Z12, who 
derived a number of important formulae for determining 
the properties of gas motions in the ICM from measure¬ 
ments of the shift and width of spectral lines. Our par¬ 
ticular case of the Coma cluster allows us to make some 
further simplifying assumptions, as we will note below. 


2.1. Cluster Model 


We construct a simple model of the Coma cluster as¬ 
suming thaWtiiS_§i®£fl 2 S_JlllSlb 6 L„^£Il®ity_Jollows a j3- 
model (jCavaliere fc Fnsco-Femiaii3ll976l Il978[ l: 


ne(r) = no 




( 1 ) 


with no = 0.003 cm“^, Tc = 300 kpc, and /3 = 2/3. We 
assume that the gas is isothermal with T = 8 keV and 
has constant metallicity Z = 0.3 Zq. These parameters 
are chosen to provide a simple model with which to per¬ 
form our calculations, and are broadly consist ent with 
/3-model fits to the Coma cluster emission (e .g., iHnghesI 
II989tlfeiel et al.lll992l : l^umann et al.ll2003ll 

Without any loss of generality, throughout this work 
we will assume the line of sight is along the 2 : direction, 
the 3D position vector x = (x, y, z), and the 3D wavevec- 
tor k = {kx, ky, kz )0 The 3D velocity field Vzi^) is as¬ 
sumed to be turbulent with an underlying power spec¬ 
trum 

PsDik) = |h.(k)p = (2) 

where Vz is the Fourier transform oi Vz, k — |k|, ko and 
ki are cutoffs of the power spectrum at high and low 


^ Here and throughout we adopt the convention that wavenum¬ 
bers k and distance scales r are related by k = 1/r. 
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Fig. 1. — Theoretical energy spectra and corresponding structure functions for a Coma-like cluster, calculated without measurement or 
“cosmic variance” uncertainty for different large-scale cutoffs (injection scales) of the velocity power spectrum. The dissipation scale for all 
curves is 20 kpc, and all spectra have a = -11/3. Colored dot-dashed lines indicate the values of the cutoff scale for each of the different 
energy spectra. The curves have been normalized by the average velocity dispersion which is the integral of over all k. Dashed 

black lines indicate the lowest and highest scales resolvable by the pointing configurations in Figure [3] 


wavenumber, respectively, Cn is a normalization con¬ 
stant, and a is the spectral index. The corresponding 
energy spectrum is given by: 

E{k) ^ P3uik)k‘^ (3) 

For a Kolmogorov power spectrum, a = —11/3 results 
in the familiar behavior for the energy spectrum 

in the inertial range of Equation [2j 


2 .2. Computing Line of Sight Velocities 


The emission-weighted line of sight velocity field is 
given by: 


Vz{x,y) 


J V2(x)e(x)dz 

fe(x)dz 


(4) 


where e is the X-ray volume emissivity neUH-l^iT, Z), 
with rie, uh the electron and proton number densities, T 
the plasma temperature, and Z the metallicity. The as¬ 
sumptions of isothermality and constant metallicity ren¬ 
der Vz dependent on neUn alone. To simplify the equa¬ 
tions, we set (as in Z12): 


resulting in 


w(x) 


e(x) 

f €(x)dz 


(5) 


Vz(x,y) 


= J Vz(x)uj{x)dz 


( 6 ) 


The line of sight velocity dispersion is similarly computed 
by: 

wl{x,y) = j v'^^{x)uj{x)dz - v‘^^{x,y) (7) 


2.3. 2D Power Spectrum and Structure Function 

Following Z12 (see their Section 3, Equa t ion 4 ; also 
Appendices A and B) and IChurazov et al.l (|2012ll (see 
their Section 3), a relationship between the power spec¬ 
trum of the 3D velocity field and that of the 2D velocity 
field can be derived: 

P2u{kx,ky;x,y) = J P^ri{\i)P^{kz-,x,y)dkz (8) 

where P 2 D is the power spectrum of the line of sight 
velocity field Vz and is the power spectrum of the 
normalized emission measure along the z-axis. They also 
pointed out that for a large wavenumbers k, the relation¬ 
ship between the 2D and 3D velocity power spectra is ap¬ 
proximately independent of the integral over the power 
spectrum of the emission measure (cf. Section 6 of Z12, 
in particular Equations 17-19 and Figure 8): 

P2D(fc) ~ P3D(fc) J Pe{kz;x,y)dkz = KPariik) (9) 

e.g., the 2D power spectrum is the same as the 3D power 
spectrum apart from a normalization constant, which is 
easily computed for our model cluster (see Appendix lAl 
for a derivation). 

In our analysis, we will use the simulated velocity fields 
to compute the second-order structure function: 

SF(r) = (|f;^(x + r)-h^(x)p), (10) 

where x ^md r are two-dimensional position vectors on 
the sky, r = |r|, and () indicates an average over pairs 
of points with the same separation r. This quantity has 
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Fig. 2.— Theoretical energy spectra and corresponding structure functions for a Coma-like cluster, calculated without measurement or 
“cosmic variance” uncertainty for different small-scale cutoffs (dissipation scales) of the velocity power spectrum. The injection scale for all 
curves is 1000 kpc, and all spectra have a = -11/3. Colored dot-dashed lines indicate the values of the cutoff scale for each of the different 
energy spectra. The curves have been normalized by the average velocity dispersion which is the integral of P 3 J 3 over all k. Dashed 

black lines indicate the lowest and highest scales resolvable by the pointing configurations in Figure [3] 


the advantage that it can be computed easily regardless 
of the spatial shape of the region under consideration. 
The second-order structure function is directly related 
to the power spectrum of the projected velocity field (see 
Appendix [B] for a derivation): 

poo 

SF{r) = Att / [1 — Jo{27rkr)] P 2 B{k)kdk (11) 

Jo 

where Jq is a Bessel function of the first kind. 

If we only consider the power-law form of the power 
spectrum, we can derive a scaling for the structure func¬ 
tion in the inertial range. Neglecting the variation in Jq, 
and assuming a characteristic scale of kc = A/r for the 
turbulent cascade, we have 

poo 

SF{r)(x / P 2 Bkdk (X k^+^ (X F, (12) 

Jkc 

with 7 = — (a + 2). For a Kolmogorov spectrum with 
a = -11/3, 7 = 5/3. 

To determine what form we may expect SF(r) to take 
for different forms of the power spectrum P 2 D{k), we will 
numerically integrate Eauation llll assuming a Coma-like 
cluster and ignoring sources of error. We choose to ex¬ 
amine two sets of models, one with different cutoffs of 
the power spectrum at small scales (fcp = l//min)) and 
another with different cutoffs of the power spectrum at 
large scales (fci = 1/Zmax)- The former represents dif¬ 
ferent possible values of the dissipation scale of the tur¬ 
bulent motions, whereas the latter represents different 
possible values of the injection scale of the turbulent mo¬ 
tions. Each set will have the same cutoff scale at the 


opposite end. Figures [T] and [5] show the energy spectra 
and structure functions calculated for the large-scale and 
small-scale cutoffs, respectively, assuming a Kolmogorov 
slope of a = —11/3. The black dashed lines in the fig¬ 
ure indicate the range of length scales sampled by the 
observations we will simulate in Section 12.51 The struc¬ 
ture functions with the small-scale cutoffs (left panel) all 
have a very similar shape over the relevant length scales, 
diverging only slightly toward the smaller scales. On the 
other hand, the curves with the large-scale cutoffs ex¬ 
hibit large differences in shape over the relevant length 
scales. As the cutoff scale (max is decreased, the struc¬ 
ture function becomes flatter, and with (max = 100 kpc 
the curve is essentially constant over all sampled length 
scales. The structure functions are normalized by the 
average velocity dispersion (which will be measured in¬ 
dependently from the line width), in order to illustrate 
the extra information contained in the shape of the struc¬ 
ture function itself. 

2.4. Defining Parameter Spaces 

To decide on our chosen values of the various input 
parameters for our velocity power spectra, we draw on 
both theoretical and observational constraints specific to 
the Coma cluster. 

We determine our range for the injection scale based 
on observati ons of the Coma c luster and theoretical con¬ 
siderations. IChurazov et al.l (|2012D used Chandra and 
XMM-Newton observations to analyze surface brightness 
fluctuations in Coma, deriving a 3D power spectrum of 
density fluctuations. On large scales, the spectrum of 
fluctuations begins to turn over for length scales larger 
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than Zmax 300 — 500 kpc (their Figures 14 and 15). If 
the density fluctuations are associated with turbulence, 
we can identify this scale with the i njection scale of the 
gas m otion. Coma’s large radio halo (jBrown fc R.udni^ 
ImID and lack of a central cool core are evidence for re¬ 
cent merging activit y, which can drive gas motions on 
scale s up to a Mpc faicker fc Sarazinll2001l : iPoole et'al] 
I2006I1 . A rece nt simulation of the formati on of a Coma¬ 
like cluster by iMiniatifc BeresnvakI (l2015f) estimated an 
injection scale of ^1 Mpc. Galaxies moving through 
the ICM can drive turbulence at smaller scales, around 
^100 kpc. From these considerations, we adopt a wide 
range of injection scales from 100 kpc up to 1000 kpc for 
our simulations. 

For the dissipation scale, we may use the iKolmogorovI 
(jl941fl formalism to estimate what a reasonable expecta¬ 
tion of its value may be. Assuming turbulence is injected 
at scale Zmax with velocity Umax, we may estimate the en¬ 
ergy injection rate as 

e ^ ^max/^max (13) 

this energy will be dissipated at the scale Zmin with a 
characteristic velocity Umin given by 

^minUmin ^ ^ ( 1 ^) 

where v is the kinematic coefficient of viscosity. Since 
in the inertial range the energy cascade rate is a con¬ 
stant, we may combine Equations [13] and [14] to give the 
dissipation scale as 



If we assume the viscosity is given by the Spitzer ([Spitzed 
1196211 formulation, given the properties of the plasma 
given in Section 12.11 and the range of velocities and in¬ 
jection scales that we assume, we derive a range for the 
dissipation scale of the turbulence of Zmin 10 — 50 kpc. 
The entire range for the dissipation scale is smaller than 
or comparable to the Astro-H resolution at the redshift 
of Coma, implying that it will be difficult to constrain 
this parameter. For this reason, we have not only cho¬ 
sen values for the dissipation scale ^min within this range, 
but also somewhat higher (100-200 kpc), to determine if 
Astro-H would be able to distinguish even extreme values 
of this scale from more physically sensible ones. 

The choice of an isotropic Spitzer viscosity constitutes 
a likely upper limit on the viscosity of the ICM. Given 
that the cluster plasma is weakly magnetized, the ion 
mean free path is much larger than it s Larmor radius , 
which renders the viscosity anisotropic (|Braginskiilll965h . 
Under the assumption of an isotropically tangled mag¬ 
netic field, thi s will suppress the average viscosity by 
a factor of 5 (jNulsen fc McNamaral [2013]) . A proper 
treatment of turbulent dissipation in such a magnetized 
medium would take into account the fact that only fast 
and slow magnetosonic modes will be damped by Bragin- 
skii viscosity, whereas Alfve nic modes will be undamped 
(jBrunetti fc LazarianI l20(T7h . The effective viscosity of 
the ICM on small scales may depend on w hether or 
not th e turbulence is super or sub-Alfvenic (jLazarianl 
l2006al[bl l. Finally, microscale plasma instabilities (such 
as firehose and mirror) may limit viscous damping 


caused by ion collisions, reduc jng the viscosity and 
the d i ssipation scale even further (iLazaria n fc Ber esnvakI 
20061: ISch ekoc hihin fc Cowle v |_[2006t iKunz et al.l 12014 
Mogavero fc SchekochihinI [^1^ These complex issues 
are beyond the scope of this paper, and beyond the ca¬ 
pabilities of Astro-H to explore, so we use the Spitzer 
viscosity as a useful benchmark. 

Finally, in choosing the spectral index a of the turbu¬ 
lent cascade, we rely on different theoretical descriptions 
of turbulence that are relevant for galaxy clusters. If the 
turbulence is primarily hydrodynamic and incompress¬ 
ible, we expect the standard Kolmogorov scaling with 
7 = 5/3 (as defined in Equation IT^ . If the turbulence is 
compressible with weak shocks in an otherwise subsonic 
velocity structure, we expect a Burgers scaling with 7 
= 2 . Alternatively, if the turbulence is magentohydro- 
dynamical in natu re, we expect a Kraichnan spe ctrum 
with 7 = 3/2 (see [Brandenburg fc Lazarianll2013L for a 
review). 

2.5. Generating Simulated Velocity Fields and their 
Structure Functions 

Using the form of the 3D power spectrum defined by 
Equation[21 we can generate realizations of the 3D veloc¬ 
ity field. Each realization of the velocity field is set up on 
a uniformly gridded rectangular domain with a projected 
area on the sky of size L^xLy = 15'x 15', or Ri 414 kpc on 
a side, and a line of sight width = 2070 kpc. The do¬ 
main is subdivided into with n^xuyXnz = 40 x 40 x 200 
cells, resulting in a cell width of r; 10 kpc on a side. This 
resolution is sufficient for our purposes, as it is several 
times smaller than the Astro-H resolution at the distance 
of Coma, approximately 40 kpc. 

A given realization of the velocity field in the z- 
direction is generated by setting up a Gaussian ran¬ 
dom field (we follow a _ pro cedure similar to that 

found in iMurgia et al.l [2004 iRuszkowski et al.l 120071 : 
iRuszkowski fc OhI 120101 : iZuHone et al.l I2011L among 
other works). We begin by computing the Eourier trans¬ 
form of the 3D velocity field, Uz(k) = ve’''^, where the 
amplitude v and the phase (/ for all grid points in fc-space 
are drawn randomly from a Rayleigh distribution: 

V ( \ dd) 

P{v, (j))dvd<p = ^ exp j dv— (16) 

where = P^D{k)/2. The normalization Cn of the 
power spectrum is set by requiring that the line of sight 
velocity dispersion {wD in the cluster averaged over a 
number of realizations is equal to a fiducial value given 
by Wzfi = -MzCs, where j\4z is the Mach number of the 
velocity dispersion, which will be varied for different re¬ 
alizations of the structure function, and Cg is the sound 
speed of the thermal gas. This value is related to the 
power spectrum by (Z12, Equation 5): 

= J P3D(k)[l - Pg{kz)]d^'k (17) 

The inverse Fourier transform of 72 (k) yields the 3D 
velocity field Vz(x). This field is then projected along 
the line of sight in the manner of equations [S] and [7] to 
produce 2D maps of the line of sight velocity (the line 
centroid) Vz(x, y) and standard deviation (the line width) 
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Fig. 3. — Configurations over which the structure function is calculated, shown as combinations of red 3’x3’ squares, overlaid on a 
sample line centroid field. Contours in the upper-left panel show the X-ray surface brightness distribution in units of ^o, the central surface 
brightness. 


Wz{x,y). These maps are then reblocked by a factor of 
4, resulting in pixels of width 1.5’ (~ 41.4 kpc), which is 
approximately the smallest angular scale that SXS will 
be able to resolve. We experimented with other options 
for reblocking the image, and found that this value rep¬ 
resents the best compromise between reducing the sta¬ 
tistical errors on the line shift and width and resolving 
power on small scales adequately. We did not include the 
effect of PSF scattering (recall for SXS, the PSF FWHM 
^1’); it will be mild for our chosen pixel size and the flat 
brightness distribution of a Coma-like cluster, and will 
be easily accounted for in the analysis of the real data. 

In our analysis, the relevant quantities related to the 
velocity field will be computed within several different 
configurations of multiple pointings of width 3’ (the SXS 
field of view). They consist of the following (seen in 
Figure[3]): a single “strip” shape 5 pointings (15’) across, 
a “big cross” shape with the pointings spaced 6’ apart, a 
“small cross” shape with the pointings spaced 3’ apart, 
a “checkerboard” conhguration that adds four pointings 
at the corners of the central pointing to the “big cross” 
configuration, and a “fill” configuration that nearly fills 
the space of the simulated data. In the “checkerboard” 
configuration, we will assume the same total exposure 
time summed across the pointings as for the “big cross” 
configuration, resulting in a reduction in exposure time 
per pointing of roughly half. For the “fill” configuration, 
we will assume the same exposure per SXS pointing as the 
“strip” and “big cross” configurations, resulting in a total 
exposure 4.2 times larger. All of these configurations fit 


within the core radius of Coma, where the density is 
still relatively fiat and the effects of stratification will be 
minor. For this reason, we have not explored pointing 
configurations with baselines larger than these. Also, 
the range of baselines these pointing configurations probe 
represents the most likely range of injection scales from 
a physical point of view (see above). 

From our maps of the line centroid, we compute the 
structure function by taking velocity differences between 
pairs of 1.5’x 1.5’ pixels and averaging these differences 
within bins at length scale r: 

SF{r) = {Sij)r = {{Vi - Vj f)r (18) 

The discrete nature of our maps results in a small num¬ 
ber of unique distance scale values between pixels, which 
suggests a simple binning scheme where similar length 
scales are grouped together. We have experimented with 
the number and position of bins and have found that our 
choice of 4-5 bins per pointing configuration, determined 
by the similarity of length scales, is adequate to resolve 
the shape of the structure function. 

2 .6. Characterization of Sources of Error 

In this section, we describe the sources of error affect¬ 
ing the measurement of the second-order structure func¬ 
tion from line shifts measured by Astro-H. In the discus¬ 
sion that follows, we only present examples where the 
injection scale of the turbulent motions is varied, since 
the variances between these different velocity fields are 

































































7 


r (arcmin) r (arcmin) r (arcmin) 

1 2 3 5 10 15 1 2 3 5 10 15 1 2 3 5 10 15 




Fig. 4.— Mean and variance of the structure functions for different large-scale cutoffs for Adz = 0.3 and the five different pointing 
configurations, computed from 100 realizations of the velocity field. No measurement errors are included, the error bars represent the 
intrinsic “cosmic variance” in the structure function. The analytic prediction is plotted with dashed lines for comparison. 


more pronounced than when the dissipation scale of the 
turbulence is varied. 

2.7. Cosmic Variance 

The structure function SF(r) (and the equivalent 
power spectrum P{k)) represents an average over all pos¬ 
sible random realizations of the underlying velocity dis¬ 
tribution function. For a given realization of a velocity 
field, a structure function may be computed from the ve¬ 
locity differences over all length scales within the field, 
but this will be only a single realization of the underlying 
distribution. A real cluster in the sky will similarly pro¬ 
duce just a single realization of the structure function. 
To represent this “cosmic variance,” for each of the dif¬ 
ferent models of the underlying power spectrum, we gen¬ 
erate 100 realizations of the velocity field. For the fc-th 
velocity field, we compute the structure function SFk{r) 
within the chosen pointing configuration, and then com¬ 
pute the mean and variance of all the SFk{r). 

Figure S] shows an example computation of the mean 
and 1-cr variance of the structure function calculated for 
three values of the injection scale /max, with all other pa¬ 
rameters held fixed. No measurement errors are included 
at this stage, only those from “cosmic variance.” This 
suggests that distinguishing between different underlying 
models will be limited by this variance. As we will see, 
the importance of the cosmic variance is different for dif¬ 
ferent choices of the coverage of velocity measurements 
from the cluster emission. Needless to say, with more 


pointings, the cosmic variance errors are reduced. This 
is especially true in the “fill” configuration. However, 
this configuration would require a large and expensive 
survery. Below, we will look for a compromise between 
the need to sample more sight lines and the need to limit 
the total required exposure. 


2 .8. Measurement Errors 

Once a given realization of the velocity field is gen¬ 
erated, the effect of statistical and systematic errors on 
the measurement of the line shift and width must also 
be taken into account. We simulate statistical errors on 
the line shift by adding a normally distributed velocity 
component with standard deviation Ustat to each pixel, 
which depends on the line width, or the Mach number of 
the turbulence. We also estimate the statistical error on 
the line width (which will be used in Section lT5l for an 
additional constraint on the velocity field parameters). 
For a given line width, the statistical error on the line 
energy scales as roughly (where Accounts is the 

number of counts in the line), and also increases for in¬ 
creasing line width, or Wzfi- Table[T]shows the estimated 
statistical errors on the line shift and width for a single 
1.5’X 1.5’ pixel for several different values of the line of 
sight turbulent Mach number Mz, assuming an exposure 
time per pointing of ^100 ks (corresponding to roughly 
125 counts in the Fe-K line per pixel for a Coma-like 









































































































Fig. 5. — The effect of bias on the computed values of the structure function, for different large-scale cutoffs for = 0.3 and the 
“strip” pointing configuration. Dashed lines indicate the analytic expectation. Solid lines indicate the computed structure function without 
measurement errors added. Open circles indicate the computed function with the bias due to measurement errors included, and filled 
circles indicate the computed function with this bias subtracted. Left panel: Bias correction for non-zero statistical error. Right panel: 
bias correction for zero statistical error (in this case, there is no bias at the smallest scale, since the velocity differences at this scale come 
from within the same pointing, and there is no systematic error between velocity measurements at this scale, see also Appendix [O) 


TABLE 1 

Statistical Errors on the Line Shift and Width 


Mz 

Wz,o (km/s) 

^stat.shift (km/s) 

•^stat,shift/'h^z,0 (%) 

<^stat.width (^^/®) 

^stat.width/‘h^z,0 (%) 

0.1 

146 

33 

22.6 

50 

34.2 

0.2 

292 

59 

20.2 

71 

24.3 

0.3 

438 

86 

19.6 

97 

19.6 

0.5 

729 

142 

19.5 

153 

21.0 


cluster) 0 Additionally, to take into account a conser¬ 
vative systematic error on the line energy from the sys¬ 
tematic uncertainty of SXS gain stability, we add a nor¬ 
mally distributed velocity component of fisys = 60 km/s 
(equivalent to ^1.3 eV at E’ ^ 6.5 keV) to the veloc¬ 
ity. The added systematic uncertainty is constant for a 
given 3’x3’ pointing, e.g., for velocity differences within 
the same pointing, the systematic error cancels out. For 
each of the 100 realizations of the velocity field, we com- 

® The results here are based on the statistical accuracy on the 
line velocity derived using ~125 line counts. Those counts can 
come from the He-line Fe line {E = 6.7 keV) alone, or from the 
He-like and H-like (6.9 keV) lines combined; the latter gives a sim¬ 
ilar formal accuracy for the same total counts. For the Coma gas 
temperature, the He-like to H-like line ratio is about 10:7, so the ex¬ 
posures that correspond to our simulated statistical accuracy would 
depend on whether the observer chooses to combine the two lines 
('^ 100 ks is for the use of the He-like line alone). This will require 
a judgment on whether the two lines sample the same gas along the 
1 . 0 .s. (e.g., there are no projected hot clouds, etc.) based on a more 
detailed study of the Coma spectrum. Here we assume that the 
line emissivity (per unit emission measure) is uniform throughout 
the cluster. 


pute 100 realizations of the measurement error, resulting 
in a total of 10^ realizations of the velocity field for each 
combination of the input parameters. These realizations 
will determine the scatter of the structure function due 
to both “cosmic variance” and “measurement” errors. 

The addition of measurement errors necessarily intro¬ 
duces a bias in the computation of the structure function. 
Since these velocity differences add in quadrature when 
we compute the structure function, the end result is that 
the observed structure function SF'[r) is biased upward 
from the “true” value SF{r). Under the assumptions of 
uniform CTstat and CTsys per pointing, we can compute the 
bias to be (see Appendix ICl) : 

SF'{r) = SF{r) + 2aLt + (19) 

where for the smallest length scales (within the 3’x3’ 
SXS field of view), a^ys = 0. 

Figure [S] shows the computed structure functions with 
and without measurement errors, to illustrate the effect 
of the bias, as well as the structure function with the bias 
subtracted, following Equation [191 First, as seen in the 
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Fig. 6. — The ratio of measurement error to cosmic variance error only for the “strip,” “big cross,” “checkerboard,” and “fill” configurations 
for structure functions with different injection scales. The “strip,” “big cross,” and “checkerboard” configurations all have the same total 
exposure, whereas the “fill” configuration has the same exposure per pointing as “strip,” for a longer total exposure. 
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Fig. 7.— Structure functions for input power spectra with different dissipation scales /mim for J^z = 0.3 and three different pointing 
configurations. Colored dot-dash lines indicate the locations of the cutoff scales. 


figure, the bias is most serious for the structure function 
with a small injection scale /max; where the true velocity 
differences are small, and dominated by the measure¬ 
ment errors. Second, the bias-subtracted curve agrees 
well with the analytic expectation, as well as the com¬ 
puted curve generated without measurement errors. We 
present two cases, one where the statistical error is set to 
the expected value for our selected exposure time, and 
another with vanishing statistical error, in the limit of 
large exposure time (so only the systematic errors, which 
are independent of exposure time, are included). We see 
from this comparison that the dominant contribution to 
the bias on the structure function is from the statistical 
error. All of our results from this point on are presented 
with this bias (both statistical and systematic) corrected. 

The measurement errors will increase the variance of 
the structure function beyond that expected from “cos¬ 
mic variance” in the velocity field. Figure [5] shows the 
ratio of measurement and cosmic variance error for differ¬ 
ent pointing configurations and different injection scales. 
The velocity field with an injection scale of 100 kpc is 
most severely affected by the measurement errors, so 


much so that they dominate the total error. For the 
other injection scales, the measurement error is typically 
comparable to the cosmic variance error, for our partic¬ 
ular choice of exposures. In particular, for the “strip,” 
“big cross,” and “checkerboard” configurations, which 
divide the same total exposure differently among several 
offsets, these errors are close for the mid-range value of 
the injection scale, which means that the exposures are 
near-optimal. 

3. RESULTS 

3.1. Small-scale Cutoffs/Dissipation Scales 

Figure [7] shows the computed structure functions for 
small-scale cutoffs of 20, 50, 100, and 200 kpc for the 
“strip,” “small cross,” and “fill” pointing configurations 
and for a value of Alz = 0.3, all assuming an exposure 
time of ^100 ks per pointing (so “fill” has 4.2 times the 
total exposure of the other two configurations). In each 
of the figures, the vertical error bars correspond to the 
1 -cr errors (including both the cosmic variance and mea¬ 
surement errors) and the horizontal error bars delineate 
the bin sizes. In all of the pointing configurations, all of 
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Fig. 8 .— Structure functions for input power spectra with different injection scales Zmax, for Jv[^ = 0.3 and four different pointing 
configurations. Colored dot-dash lines indicate the locations of the cutoff scales. 


these curves are indistinguishable within the 1-cr errors. 
In addition, the scales below ^40 kpc are not resolved 
by the 1.5’ pixels. Unless the dissipation scale is much 
larger than 200 kpc (an unexpected scenario, given the 
theoretical motivations for the expected range of this pa¬ 
rameter; see Section [5^, Astro-Hw\\\ not be able to dis¬ 
tinguish the corresponding structure function from one 
with a smaller dissipation scale. 

3.2. Large-scale Cutoffs/Injection Scales 

Figure [5] shows the computed structure functions for 
large-scale cutoffs of 100, 200, 300, 500, and 1000 kpc for 
the different pointing configurations, for a value of Mz 
= 0.3, assuming a ~500 ks total exposure for each con¬ 


figuration, with the exception of the “fill” configuration, 
which has a total exposure of ^2100 ks. This shows that 
the possibility exists for distinguishing between structure 
functions with different injection scales, due to the very 
different shapes of the curves. The sensitivity to the 
injection scale depends on the scale of interest and the 
configuration of the pointings. 

Distinguishing between the curves at small separations 
{r/ 100 kpc) is difficult due to the proximity of the 
curves, regardless of the pointing configuration. The 
curves will be most easily distinguished at large spatial 
scales r. At these scales, the “big cross” configuration 
(top-right panel) allows for the differences in the struc¬ 
ture functions to be better constrained than the “strip” 
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configuration, due to the increased number of pairs of 
points at large distances. On the other hand, the “strip” 
pointing allows for the shape of the structure function 
in the intermediate range of r (^80-150 kpc, 3-5’) to 
be better constrained, due to the availability of distance 
measurements within this range for this pointing. This 
is the range of length scales where the structure function 
from velocity fields with = 200 and 300 kpc tran¬ 
sitions from a power-law dependence on r to a constant 
value. 

The “fill” configuration (lower-right panel) combines 
both of the strengths of these two pointings, and im¬ 
proves the resolving power at all scales, making it much 
easier to distinguish between the curves. However, it re¬ 
quires 4.2 times the exposure of the other configurations. 
The “checkerboard” configuration represents a compro¬ 
mise between the “big cross” and “fill” configurations, by 
adding four pointings at the corners of the central point¬ 
ing to increase the number of measurements at small- 
to-intermediate scales. Here, we have assumed that the 
total exposure time (summed over the pointings) for the 
“checkerboard” is the same as the total for “big cross” 
(^500 ks), so that each individual 3’x3’ pointing has 
an exposure time of ^55 ks (corresponding to ^70 line 
counts per 1.5’x 1.5’ pixel). From the perspective of be¬ 
ing able to distinguish between the structure function 
curves at intermediate to large length scales, “checker¬ 
board” is a modest improvement over that of “big cross” 
(see also Section EH). We note that we assume the sta¬ 
tistical error on the line shift scales as \/fVcounts, at even 
low values of Accounts) this will have to be verified using 
the early (real) SXS data. 

For smaller injection scales, the inertial range of the 
structure function will be difficult to discern due to the 
inability to resolve length scales below ^40 kpc, the 
smallest length scale we resolve. The ^max = 100 kpc 
and 200 kpc curves are essentially fiat over the entire 
range of r that is discernable by any of the chosen point¬ 
ing configurations. Additionally, these curves are most 
affected by the measurement errors, due to the smaller 
velocity differences. In particular, the Zmax = 100 kpc 
curve is often nearly consistent with zero due to the un¬ 
certainty from the combined effect of the statistical and 
systematic errors on the line shift. 

3.3. Different Power-law Slopes 

It is of interest to determine the slope of the structure 
function a in the inertial range (provided this range can 
be discerned; that is, there is enough of a dynamic range 
between the injection and dissipation scales), since this 
is directly related to the power-law slope of the under¬ 
lying velocity power spectrum, the value of which de¬ 
pends on the physics of the turbulent medium. Figure 
El shows a comparison of structure functions for three 
physica lly m otivated values for the slope a discussed in 
Section 12.41 with the cutoff values of /min = 20 kpc and 
/max = 1000 kpc held fixed for the different input power 
spectra, assuming ^100 ks exposure per pointing in the 
“fill” configuration. Similarly to the situation with the 
small-scale cutoffs, the curves are too similar to be dis¬ 
tinguished within the errors. Despite the qualitatively 
different physics that produces the different power spec¬ 
tra, the associated range of spectral indices is not wide 
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Fig. 9.— Structure functions for input power spectra with differ¬ 
ent power-law slopes a, for J\A.z = 0.3, and the “fill” pointing con¬ 
figuration. For all the curves, Imin = 20 kpc and /max = 1000 kpc. 

enough to be discerned by this technique. 

3.4. Different Mach Numbers 

The broadening of emission lines will provide a measure 
of the strength of the turbulence in the velocity compo¬ 
nent along the line of sight. An increase in the velocity 
scale of the turbulence will correspond to an increase in 
the overall normalization of the power spectrum as well 
as the measured structure function. Stronger turbulence 
also implies higher velocities, which will make it easier 
to measure velocity differences. 

On the other hand, for a constant number of counts 
in the line, the statistical error on the line shift will be 
roughly proportional to the measured width of the emis¬ 
sion line, which contains contributions from thermal, in¬ 
strumental, and turbulent broadening: 

^obs = + ■a’therm + ^turb (20) 

In our situation, the first two terms will be roughly con¬ 
stant, implying for low Wturb the error on the line shift 
will be approximately constant, and increase for stronger 
turbulence. The systematic error on the line shift is in¬ 
strumental in nature and will be independent of the tur¬ 
bulent velocity strength. 

Figure (Tn] shows the observed structure functions for 
different values of the line of sight Mach number Aiz, 
assuming ^100 ks exposure per pointing. As the over¬ 
all strength of the turbulence decreases, it will be more 
difficult to distinguish velocity fields with different injec¬ 
tion scales from each other, due to the fact that as the 
magnitude of the velocity difference decreases, the mea¬ 
surement errors on the line shift are roughly constant. 
For turbulent Mach numbers A4z < 0.2 in Coma, we are 
unlikely to be able to constrain the injection scale of the 
turbulence. 

3.5. Estimating the Injection Scale of Turbulence 
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Fig. 10.— Structure functions for input power spectra with different injection scales /max? for different Mach numbers A ^2 and the “strip” 
configuration. Colored dot-dash lines indicate the locations of the cutoff scales. 


Our results indicate that a measurement of the 2D ve¬ 
locity structure function of the Coma cluster by Astro-H 
may be useful for determining the injection scale of the 
turbulent motions. To determine more precisely the con¬ 
straints that future observations can put on this param¬ 
eter, we jointly fit the model for the structure function 
(Equation If II) and the average velocity dispersion within 
the total configuration (Equation ITTl) for the parameters 
Cn (the power spectrum normalization) and (max (the in¬ 
jection scale) to the mean structure function and average 
velocity dispersion computed for each pointing configu¬ 
ration. Note that the velocity dispersion, which we have 
used to define the Mach number provides an addi¬ 
tional constraint on the parameters Cn and (max, but is 
still dependent on both (see Equation ITTl) . 

We perform the fit by setting up a finely spaced 2D grid 
in Cn and (max, and for each grid point we calculate the 
statistic. We compute the 1-cr errors on the structure 
function and average velocity dispersion from our sample 
of realizations of each, which include both measurement 
and “cosmic variance” errors. All other parameter values 
are held fixed, at a = —11/3 and (min = 20 kpc (other 
physically sensible values for these parameters will not 
change our conclusions, given our insensitivity to these 
parameters, see Sections l3 .1 1 and 183)) . Under these con¬ 
straints, for a measured value of the average velocity dis¬ 
persion {wz) {fx Aiz), a given value of either Cn or (max 
uniquely determines the value of the other parameter. 

Eigure [TT] shows plots of the 1-parameter 68% (y^ < 
Xmin + 1) confidence regions of the parameters (max and 
Cn, for .Adz = 0.3 and different values of the true (max, for 
the different pointing configurations. Table [2] shows the 
best-fit (max and l-tr errors for each of the configurations. 

In general, we find that there is a parameter degener¬ 
acy between the normalization Cn and the injection scale 
(max for a given A4z (see the long banana-shaped confi¬ 
dence regions in Figure El). Figure El illustrates this, 
by showing an example best-fit curve and two curves at 
extreme values of the parameters Cn and (max, at the ex¬ 
treme ends of the two-parameter 68% confidence region 
of the data (Xmin + 2-3), both of which fit the structure 
function and the velocity dispersion. This figure illus- 


TABLE 2 

Constraints on /max 


/max (kpc) 


Estimated /max (kpc) 



strip 

big cross 

fill 

checkerboard 

100 

lObtlt 

103/1® 

103/18 

103/^1 

200 

209t(?i* 

2051“® 

205/4° 

201 /8° 

300 

329ifj^ 

329/f?f 

322/g 

32214844 

500 

623™ 

582±®13 

552/241 

562/442 

1000 

1339+A|g’ 

1208™ 

1062/|8| 

1120115180 


trates clearly how a lower Cn is compensated by a larger 
(max, and vice-versa. 

The “strip” configuration (upper-left panel) does the 
poorest job of constraining these parameters, as shown 
by the overlapping confidence regions in Figure ITT] The 
“big cross” configuration (upper-right panel) fares bet¬ 
ter, due to the larger number of baselines at large dis¬ 
tance scales, and provides tighter limits on (max- As seen 
in Figure El and Table [21 constraining (max at values 
larger than those that are covered by our region of inter¬ 
est (^400 kpc wide) is difficult. Such cluster-scale eddies 
might be caused by a recent major merger. 

It is readily apparent from Figure El (and not unex¬ 
pected) that the “fill” configuration with the same expo¬ 
sure time per pointing (and much longer total exposure, 
lower-right panel) does the best job of constraining the 
injection scale and normalization of the turbulent cas¬ 
cade. However, this is a very expensive survey. The 
“checkerboard” configuration (lower-left panel), with its 
lower exposure time per pointing but more pointings, 
offers a modest improvement over “big cross” in terms 
of constraining the injection scale. However, it relies 
on measuring velocities based on a small number of line 
counts, the possibility of which needs to be demonstrated 
with real data. 

3.6. Multiple Injection Scales 

In reality, the distribution of turbulent velocities in 
the Coma cluster is likely to be more complicated than 
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Fig. 11.— 1-parameter 68% confidence limits on the Zmax and Cn parameters for four different pointing configurations and different 
values of Imax- Crosses mark the positions of the true parameter values. 


the simple power-law models we have investigated in this 
work. One likely possibility, motivated by cosmological 
simulations of merging clusters, is that there will be more 
than one scale at which significant driving of turublent 
motions is occurring. In the Coma cluster, there is an 
ongoing cluster merger which is driving turbulence at 
large scales (a few hundred kpc), whereas the motion of 
galaxies within the ICM will stir the gas at scales of 100 
kpc or less. 

Using numerical simulations, lYoo fc Choi (larfli ex¬ 
plored the effect of driving turbulence on multiple length 
scales. They found in the case of driving on two scales 
that the relative height of the peaks of the power spec¬ 
tra and the separation of the driving scales is crucial to 


whether or not the two components can be distinguished. 
In particular, it will be easiest to discern the presence of 
driving on two length scales if the two peaks are roughly 
equal in height and the driving scales are well-separated. 

To investigate the ability of Astro-H to distinguish mul¬ 
tiple scales, we modeled an input power spectrum that 
is the sum of two single-injection scale models (shown in 
the left panel of Figure flSl). one with M — O.b and an 
injection scale of Zmax = 100 kpc and another with A4 = 
0.2 and an injection scale of Zmax = 1000 kpc. The green 
curve shows the total model, and the red and blue curves 
show the individual subcomponents of this model. What 
results is an energy spectrum with two resolvable peaks 
of approximately the same height, separated enough in 
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Fig. 12.— Best-fit and example curves with extreme parame¬ 
ter values (taken from the extreme ends of the 2-parameter 68% 
confidence regions) compared to the data for the “big cross” config¬ 
uration, A4z = 0.3, and /max = 300 kpc. The main panel shows the 
structure function, while the inset shows the velocity dispersion. 


length scale so that they are distinct. We use this in¬ 
put power spectrum to generate velocities in the same 
manner as our previously described models. 

The right panel of Figure [13] shows the correspond¬ 
ing structure function for this model. The dashed green 
line shows the expected curve. The structure function 
has a flattening at intermediate r that corresponds to 
the smaller injection scale, then an increase followed by 
another flattening corresponding to the larger injection 
scale. The solid lines with error bars show the computed 
structure functions from the “fill” pointing configuration, 
chosen to maximize the probability that the two injec¬ 
tion scales can be discerned. However, even in this case, 
due to the finite resolution the computed structure func¬ 
tion is unable to resolve the transition region between 
the two subcomponents at an angular scale of approx¬ 
imately 3.5’ (length scale of 100 kpc) between the two 
injection scales, and instead the behavior of the com¬ 
bined model appears indistinguishable from a power-law 
within scales 40 kpc ^ r < 300 kpc. The slope appears 
shallower than those in the single-injection models, since 
the reduction in the structure function with decreasing 
r from the /max = 1000 kpc component is partially offset 
by the contribution from the /max = 100 kpc component. 

This demonstrates that distinguishing between mul¬ 
tiple injection scales will be difficult, not only due to 
the cosmic variance and measurement uncertainty, but 
also from the limited range of length scales that Astro-H 
will be able to sample. It is also clear from Figure |T3| 
that a mission with a higher spatial resolution, such as 
Athena or Smart-X, would present a better opportunity 
for discerning the two driving scales, since the slope of 
the lower-injection scale model steepens at angular scales 
less than 1.5’ (length scales < 40 kpc). 


4. SUMMARY 

We have performed simulations of turbulent velocity 
fields in the Coma cluster, and constructed the 2D struc¬ 
ture function of the line of sight velocity field (line shift). 
These simulations bear directly on future observations 
of the Coma cluster by Astro-H. We have simulated a 
number of realistic pointing configurations and models 
for the underlying velocity field. Our main results are as 
follows: 

• It will not be feasible to distinguish between struc¬ 
ture functions with different dissipation scales, re¬ 
gardless of the number and layout of the pointings 
used, unless this scale is implausibly large (on the 
order of 100s of kpc). 

• It will be possible to measure the turbulence injec¬ 
tion scale, for any of the mapping configurations 
examined in this work, though some configurations 
resolve some length scales better than others. How¬ 
ever, this is likely only feasible for turbulence where 
the Mach number of the line of sight velocity dis¬ 
persion is Af z ^ 0.2, unless the systematic uncer¬ 
tainties are much better than anticipated. 

• It will not be feasible to distinguish between dif¬ 
ferent power-law slopes in the inertial range of the 
power spectrum, unless the slopes are significantly 
steeper or shallower than suggested by physically 
motivated models of turbulence. 

• Since the structure functions for different injection 
scales are best-distinguished by their behavior at 
large distance scales, the “big cross” pointing con¬ 
figuration is preferable to the “strip” pointing con¬ 
figuration, as it samples larger length scales better 
at the expense of intermediate scales. However, 
we find it is advantageous to split the same total 
exposure among a larger number of offsets. The 
“checkerboard” configuration places tighter con¬ 
straints on the injection scale, though these im¬ 
provements are modest. We therefore recommend 
that Coma should be mapped with either the “big 
cross” pointing at ^100 ks per pointing, or with 
the additional four pointings required to make the 
“checkerboard” configuration at ^55 ks per point¬ 
ing. The latter configuration will reduce the cos¬ 
mic variance error on the structure function, but it 
requires conhdence (beyond the formal statistical 
sense) in velocity measurements derived from such 
a low number of counts. 

Other statistics constructed from the measured veloc¬ 
ity field of Coma may also be of use in determining the 
properties of the velocity field. Z12 suggested using the 
ratio of the root mean square of the line shift to the 
velocity dispersion as a function of length scale to con¬ 
strain the injection scale of turbulence, a measure similar 
to our joint fit of the structur e fun ction and the average 
velocity dispersion in Section 13.51 We have also not ex¬ 
amined the spatial dependence of the velocity dispersion 
(line broadening), which may place additonal constraints 
on the turbulent power spectrum. 

Furthermore, the non-Gaussian shape of the turbu- 
lently broadened line profiles also encodes information 
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Fig. 13.— Driving turbulence on more than one scale. Left panel: Energy spectra for a model where the total spectrum is a sum of two 
components with different driving scales. Right panel: Expected and computed structure function for the two-scale model. Red and blue 
curves show the individual model components, whereas the green curve shows the total model. 


on th e turbulent power spectrum (jlnogamov fc: SunvaevI 
120031) , but to access this information, much longer expo¬ 
sures will be required. For well-resolved clusters such 
as Coma, we can constrain the power spectrum with 
cheaper, spatially resolved measurements, as shown in 
this paper. However, for those clusters where the Astro- 
H angular resolution will be in, the information in the 
shape of the line profile may be essential. 

Finally, though similar analyses may be performed 
with other clusters, it will be more difficult. We have 
chosen Coma because its large core radius and lack of 
a central AGN make it a relatively clean setup for this 
kind of analysis. Other nearby systems, such as Virgo 
and Perseus, have the requisite spatial resolution, but 
have steep entropy gradients and a large dynamic range 
in density and pressure within the core region. Addi¬ 
tionally, they are both sources of AGN activity and pos¬ 
sess sloshing cold fronts, which will limit the locations 
at which we can expect to get a measurement of the 
velocity field that is not contaminated by these effects. 
Perhaps most importantly, the simplifying assumptions 
of homogeneous and isotropic turbulence made in this 
work will have to be relaxed, since the significant strat- 
ifcation of the cluster atmosphere in these systems will 
prevent large eddies from developing in the radial direc¬ 
tion, whil e permitting them to dev e lop in the tangential 
direction (jZhuravleva et al1l2014al I2015D . Z12 showed 
that the velocity dispersion as a function of projected ra¬ 
dius from the cluster center is effectively a measurement 
of the structure function, though we must add a caveat 
that for clusters with highly stratified atmospheres the 
differences in scales probed by this method are signifi¬ 
cantly affected by differences in the power spectrum itself 


as a function of radius. 

Our results have shown that Astro-H will be able, for 
the first time, to directly measure the Mach number and 
the injection scale of the turbulent velocity field of a 
galaxy cluster, in this case the nearby Coma cluster. 
These are important quantities, as they (a) determine 
the diffusion rate for everything advected by the gas, 
including cosmic rays, metals, etc., and (b) the rate of 
energy flow down the turbulent cascade, which ends up 
dissipated into heat and nonthermal components. Re¬ 
solving other aspects of ICM turbulence that have bear¬ 
ing on the underlying physics will require a mission with 
finer angular resolution and larger effective area, such as 
the upcoming Athen^ or the proposed Smart-^^ These 
include resolving the turbulent dissipation scale and con¬ 
straining the spectral index of the turbulent power spec¬ 
trum, both of which depend on the microphysics of the 
plasma and the driving mechanism of the gas motions. 
The lessons learned in this work can be applied to future 
investigations assuming the characteristics of any future 
X-ray instrument. 
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APPENDIX 

A. CALCULATION OF THE NORMALIZATION OF THE 2D POWER SPECTRUM 

Section 12.31 indicated the relationship of the velocity power spectrum in two and three dimensions can be approxi¬ 
mated by a normalization constant given by 


K = 


Pe{k^;x,y)dkz 


(Al) 


where Peix,y', kz) is the power spectrum of the normalized emission a;(x) for a given position {x,y) on the sky. Our 
choice of /3 = 2/3 in Equation[T]simplifies the required integrations considerably. Setting = x^+y"^ and 
we can rewrite Equation [S] as 


w(x) 




1 + 



1 





2c3 1 

TT (c2 + 2;2)2 


(A2) 
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the Fourier transform of oj is therefore given by 

2c3 


u}{k^-x,y} = - / 


— 2T^ik, 


^^dz = e-2’^"l'=^l(27rc|/c^| + l) 


/_oo (c^ + z'^Y 

since w is even the imaginary part of the Fourier transform is zero. Finally, setting P^{x,y; kz) = |wp, we have: 

5 


/ OO pOQ 

PYkz;x,y)dkz= / e-4""l'=^l(27rc|A:,| + l)2dA:, = 

-OO j —OO 


47r-^i?2 + 


(A3) 


(A4) 


The average value of K over our region of interest (a radius of approximately 200 kpc) is ^1 Mpc 


B. DERIVATION OF THE STRUCTURE FUNCTION 
The structure function SF{r) (Equation [10]) can be rewritten as: 

SF{r) = {vl{x + r)) + {YHx)) - ‘^•{vzix + (Bl) 

Via the Weiner-Khinchin theorem, the inverse Fourier transform of the 2D power spectrum P 2 D(k) is the velocity 
autocorrelation: 

F(r) = {vz{x + r)Tz(x)) = J P 2 iy{k)Yk (B2) 

Since the velocity field is stationary, the first two terms of Equation |BT] are: 

{vlix + r)) = (vlix)) = r(0) = J P 2 B{k)(fk, (B3) 

yielding 

SF{r) = 2(r(0) - F(r)) = 2 J [l - P 2 r,{k)d^k = 2 J [1 - cos(27rk • r)] P 2 Bik)d'^k. (B4) 

where the last eq uality follows since the integral is even. If we take the integration in polar coordinates and use 
Parseval’s integral CArflcen fc: Webedl2005h . we finally have 

poo 

SF{r)=FK / [1 — Jo{2'Kkr)\P2B{k)kdk (B5) 

Jo 

where Jq is a Bessel function of the first kind of order zero. 


C. CALCULATION OF THE BIAS ON THE STRUCTURE FUNCTION DUE TO MEASUREMENT ERRORS 

For a single pair of velocities Vi, vj separated by a distance we define 

Sij = iv^-VjY (Cl) 

The structure function SF{r) is an average over the at a given length scale r. The observed velocities will include 
measurement errors, which can be modeled as additions to the pointwise velocities, giving the observed structure 
function value as 

s'ij = [{vi + 5vi) - {vj + dvj)Y (C2) 

where the Svi, Svj are normally distributed with zero mean and variances tr^, ct|. Expanding and rearranging terms, 
we have 

fiC = {vi — VjY + dvY + SvY — ^SviSvj + 25vi{vi — Vj) — 2Svj{vi — Vj) (C3) 

The total observed structure function is an average over the for a given length scale r, yielding 

YijY ~ ) Y F {dVi Y F {dVj Y (2Svi6vYr F I^j))r Ydvj{yi (C4) 

~ ((Tz Xj'} )r F {dVi Y F {dVj Y ‘^{dViY{dvYr F 2,(^dVi)‘^{dVjYiYi Tj))r 

where the second line follows from the statistical independence of the measurement errors from different locations and 
from the value of the velocity. Since {dvi)r = {dvj)r = 0, this reduces to 

{s'iYr = {{Vi - VjY)r F {dvY)r + {dVjYr (C5) 

= {Sij)r F (j'f F cr| 

Finally, since ai = aj, and assuming at a given point i af = F the sum in quadrature of the statistical and 
systematic errors, we have 

Yij)r = {sij)r F F 2crgyg (C6) 







18 


For this work, we have assumed that velocity differences on the smallest length scales (less than « 1.5’) have zero 
systematic error, to simplify the calculation of the bias. Though this is only strictly true for velocity differences on the 
same pointing, we find that including systematic errors on velocity differences at this small scale has little effect on our 
results, due to the fact that the number of differences within pointings is much larger than the number of differences 
across pointings for any of our setups. 



